Towards linear phononics and nonlocality tests in ion traps 
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We explore the possibility to manipulate 'massive', i.e. motional, degrees of freedom of trapped ions. In 
particular, we demonstrate that, if local control of the trapping frequencies is achieved, one can reproduce 
the full toolbox of linear optics on radial modes. Furthermore, assuming only global control of the trapping 
potential, we show that unprecedented degrees of continuous variable entanglement can be obtained and that 
nonlocality tests with massive degrees of freedom can be carried out. 
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The last decade saw a boom in the development of exper- 
imental capabilities available for quantum information pro- 
cessing. The ability to manipulate the information of discrete 
variables encoded in polarisation, spin and atomic degrees of 
freedom has by now reached very high standards. On the 
other hand, the control of continuous variable degrees of free- 
dom is still almost exclusive to light fields in quantum optics. 
Even though quantum optical systems rely on well established 
tools and are very promising for, communication tasks, they 
also suffer from significant drawbacks. Notably, the entan- 
glement generation in such systems is strongly limited by the 
efficiency of parametric processes in nonlinear crystals; more- 
over, 'static' optical degrees of freedom - i.e. light resonating 
in cavities - are seriously affected by losses and decoherence 
over the typical dynamical time scales. 

Here, we discuss in detail the possibility of controlling the 
radial motion of trapped ions, described by continuous vari- 
able quantum degrees of freedom which we will refer to as 
radial modes 01 E3]. We highlight the remarkable potential of 
such modes in view of the refined technology that has been de- 
veloped in ion traps. In this experimental setting, we demon- 
strate that any linear optical operation can be obtained for ra- 
dial modes of trapped ions by controlling the individual ra- 
dial trapping frequencies. Furthermore we show that, even if 
only global control of the trapping potential is possible, such 
systems would outperform optical modes in both achievable 
degrees of entanglement and decoherence rates. Finally, as 
an application, we consider the violation of non locality tests 
with radial modes, and show it to be achievable with cur- 
rent technology, demonstrating the potential of such setups 
not only for information processing but also as probes of fun- 
damental physics. 

The trap - We shall consider the radial modes of n ions of 
mass m and charge ze in a linear Paul trap 0]. Let Xj and 
Pj be the position and momentum operators associated to the 
radial degree of freedom of the j-th ion, which is trapped in 
the radial direction with angular frequency ojj . In the follow- 
ing, the longitudinal trapping frequency v t will be the unit 
of frequency and will set the unit of length as well (equal 
to \J z 2 e 2 / (Aireornvf), where So 1S the dielectric constant); 
also, we shall set ft = 1. The Coulomb interaction affects 
the local radial oscillation frequencies: for convenience, let 
us then define the 'effective' local radial frequencies Vj = 

— ui\ 3 , {uj} being the equilibrium po- 



sitions of the ions in the length unit set by the longitudinal 
frequency J3l. Rescaling the canonical operators according to 
Xj = ^/mvjXj, Pj = Pj/^mUj, and grouping them together 

in a vector of operators R = (x%, . . . ,x n ,pj, . . . ,p n ) T , al- 
lows one to express the global Hamiltonian of the system in 
the harmonic approximation as the following quadratic form 



H = -R T 
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where v is a diagonal matrix: v = diag (y\, . . . , v n ), while 
the potential matrix k has diagonal entries Kjj = Vj and off- 
diagonal entries Kjk — ^/ {^/ v j v k\uj — Uk\ 3 ) for j ^ k. Let 
us also recall that the canonical commutation relations can be 
expressed as [Rj,Rk] = tiljk, where the 2n x 2n matrix 
51 has entries Vtj^ = 5 n ,k-j — $n,j-k f° r 1 < j, k < 2n, 
and that, in quantum optics, Gaussian states are defined as 
states with Gaussian characteristic function: a Gaussian state 
q is completely determined by the "covariance matrix" (CM) 
cr, with entries a jk = Tr [{Rj,R k }g]/2 - Tr [RjQ]Tr [R k g] 
and by the vector of first moments R, with components Rj = 
Tr [Rjg], in terms of the vector of canonical operators R ^ 
0. 

Linear phononics - In the first part of the paper, we shall 
assume that the trapping frequencies {ojj}, and thus {vj}, 
can be controlled locally and changed suddenly. This may be 
achieved by building small, local radial electrodes, by adding 
local optical standing waves J3], or in Penning trap arrays 
iH @]. Our first aim here is to show how, in principle, this 
control allows one to perform any arbitrary 'linear optical' 
operation on the radial modes of the ions, that is any uni- 
tary operation under which, in the Heisenberg picture, the 
vector of operators R transforms linearly: R i— > SR. The 
matrix S has to be 'symplectic', i.e. S T flS = O, to preserve 
the canonical commutation relations. Any symplectic oper- 
ation S on a system of many canonical degrees of freedom 
("modes") can be decomposed into a combination of generic 
single-mode symplectic transformations and two-mode rota- 
tions ("beam splitters", in the quantum optical terminology) 
fiol [Till . It is therefore sufficient for us to establish the pos- 
sibility of performing these subclasses of operations on our 
system of n ions by manipulating the local frequencies. Sin- 
gle qubit operations - In what follows we assume that the 
original frequencies of the ions are different but commensu- 
rate, as given by, say, Vj = jv, and that v is large enough 
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so that interaction between ions suppressed IU2I1 . Let us then 
consider the reaction of the system if the frequency of the j-th 
ion changes suddenly from Vj to atjVj, for some real ctj. The 
Heisenberg equation of motion for stj and pj can be imme- 
diately integrated in such a case, resulting into a symplectic 
transformation Sj (t) 








(2) 



with c = cos(i/jOj£) and ,s = sin(zA,a,i). The first and last 
factor of this decomposition are 'squeezing' operations in the 
quantum optical terminology, whereas the second factor is 
known as a 'phase shift' (i.e., a rotation in the single-mode 
phase space). Combinations of squeezings and phase-shifts 
make up any possible single-mode symplectic operation: we 
thus need to show that such operations can be implemented 
individually on any ion of the system in a controllable man- 
ner. 

Phase-shift - To realise a phase-shift operation on the fc-th 
ion, it is sufficient to change the frequencies of all the other 
ions in the same way, such that = 1 and aj = a ^ 1 for 
j k. As apparent from Eq. (01, after a time t a = 2-k l(yoi) 
one has Sj = I2 for j ^ k (let us recall that Vj = jv by 
assumption), whereas the oscillation of the fc-th ion will have 
acquired a phase ifk = 2irk/a (with no squeezing, as at is 
kept equal to 1). If the frequencies are switched back to the 
original values after a time t a , the net effect of the evolution 
is then analogous to an 'optical' phase-shift on the ion k. 

Squeezing - In order to squeeze the state of ion k, one can 
conversely change only the pertinent frequency, so that ^ 
1 and a.j = 1 for j ^ k. Then, after a time period t = 
2ir j v, all the other ions will have returned to the initial state, 
while ion k will be squeezed and phase-shifted according to 
Eq. (ff]). Notice that the phase-shift can always be corrected 
by applying the strategy described above. Let us remark that 
the degree of squeezing achieved in Eq. (0 depends crucially 
on the phase-shift operation, as the two squeezing operations 
act along orthogonal directions and are the inverse of each 
other. In the case = (t + h)/k for ft, G N, the phase- 
shift can be balanced by a counter-rotation of it /4 in phase 
space and the final squee zing operation is a diagonal matrix 
given by diag(afc, a;^ 1 ) 11311 . Also notice that, by placing 
the ions inside cavities, the squeezing of the massive degrees 
of freedom could be transferred to light, so that radial modes 
could act as an effective source of squeezing (and potentially 
even entanglement) for optical systems as well. 
Beam-splitters - Let us now turn to 'beam-splitting' opera- 
tions(eq|2]i between any two radial modes. To this aim, it is 
sufficient to bring two modes (hereafter labeled by j and k) 
to the same frequency v = Vj = Vk, so that the Coulomb 
interaction between them is no longer suppressed. Switch- 
ing to an interaction picture, one has the following interac- 
tion Hamiltonian between the two modes: KjkXj(t)xk(t) = 



{a.j 



t Jivt 



)(a fc e" 



e*"'), where the ladder 



operators are defined as Xj ~ (a,j + aj). If the frequency v 



is sufficiently large the rotating wave approximation applies 
to yield Kjk(a,ja\, + a^ak) ■ This Hamiltonian realises exactly 
the desired beam splitter-like evolution, resulting into a sym- 
plectic transformation which mixes Xj with Xk and Pj with pk 
(rotating such pairs equally, by the angle /%-fct). For instance, a 
'50 : 50' beam splitter is achieved after a time t = n/(AKjk). 
Since the interaction requires a change of the local frequen- 
cies it includes automatically in it a local operation, which 
may however be corrected before or after the 'beam-splitting' 
procedure. 

Summing up, we have shown that any symplectic (i.e. "lin- 
ear optical") operation, including squeezing, can be imple- 
mented for radial modes of trapped ions by a proper tuning of 
the frequencies of the microtraps. Displacement operations on 
individual ions, which shift the operators Rj by a real number, 
can also be implemented in microtraps by shifting the radial 
equilibrium position of the ion, or as in 11511 . Because the free 
evolution rotates the state of the radial modes in phase space 
(see Eq. (f2]i for ctj = 1), if the operation is carried out at the 
proper time such a shift can be implemented in any direction 
of phase space and not only in the positions Xj . The unitary 
operator displacing the canonical operators of mode j by, re- 
spectively, Xj and pj will be denoted by Dj(xj,Pj). 

These findings show that all the developments based on 
Gaussian states in the quantum optical scenario, in particu- 
lar concerning entanglement manipulation and informa- 
tion protocols Hill , could be carried over to radial modes of 
ion traps if local control is achieved. In fact, linear optical op- 
erations, complemented by displacements, correspond to all 
the unitary transformations that preserve the Gaussian charac- 
ter of the initial state. 

Notably, even non Gaussian states can be engineered in 
this setup with relative ease, either by entering the nonlinear 
regime of the Coulomb interactions or by exploiting the inter- 
nal degrees of freedom of the ions. The latter also allows for 
Gaussian and non Gaussian measurements on individual ions: 
the tomography of trapped ions, corresponding to homodyn- 
ing, was proposed in ll5Ul7lll8ll and partially realized in 11911 . 
while local number states and parity could be measured using 
the scheme suggested and realized for cavity QED in I20I1 . 
Quite remarkably, such a scheme would allow one to measure 
parity on a single copy of the state and run of the apparatus. 
In the remainder of the paper, in order to demonstrate the po- 
tential of Gaussian states of radial modes in experimentally 
accessible settings, we will consider entanglement generation 
and nonlocality tests requiring only global control of the trap- 
ping potential. 

Entanglement generation - The specific Gaussian situation 
we shall address starts off from the ground state g g of Hamil- 
tonian (Q]) - with all frequencies being equal, i.e. u>j = uji for 
1 < j < n - as the initial state (which can be well approxi- 
mated in the laboratory by cooling the system to its ground 
state 12110 . Next, the frequency is changed to ojf, so that 
the state g g will not be stationary anymore under the mod- 
ified Hamiltonian. For large uii, the initial state g contains 
very little entanglement but entanglement builds up during 
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FIG. 1: Entanglement between first and last ions of the chain as a 
function of time for initial frequency uji = 100MHz, final frequency 
uif = 2MHz, f t = 1MHz, temperature T = 294° K (correspond- 
ing to TV = 2 x 10 7 thermal phonons) and different couplings to the 
environment 7; thicker curves refer to ions 1 and 2 for n = 2, while 
thinner curves refer to ions 1 and 3 for n = 3. The curves for 7 = 
and 7 = 10~ 6 Hz are very close: decoherence is almost negligible 



for such heating rates. 



the subsequent evolution (see |22|| for an analogous scheme 
in chains of nanomechanical oscillators). Entanglement may 
be quantified by the logarithmic negativity £V = log 2 ||f?||i, 
where ||g||i stands for the trace norm of the 'partially trans- 
posed' density matrix of the considered system (with this def- 
inition, a Bell pair has E_\r = 1), which is computable for 
Gaussian states [23J. The ground state of Hamiltonian (fl~|i 
is just a Gaussian state with a block diagonal CM cr g = 
{a x © cr- 1 )/2, where a x = v 1 I 2 {v 1 / 2 kv 1 / 2 )~ 1 I 2 v 1 ' 2 , and 
vanishing first moments. After the change of potential, re- 
sulting into the new quadratic Hamiltonian H' = R T H'R 
(where the 'Hamiltonian matrix' H' is implicitly defined), 
the evolved state after a time t is a Gaussian state with CM 
given by at = St&gSj, for St = exp(Q.H't), Further- 
more, we have taken into account decoherence in an envi- 
ronment of phonons with temperature T and Toss rate' 7, 
as the master equation under such conditions admits Gaus- 
sian solutions with CM e -7t er t + (1 — e~ 7t )St<r 00 Sj , where 
a-oo = diag(i + Nj, \ + Nj) and Nj = l/(c^/ T - 1) is the 
number of thermal phonons at frequency Vj (setting kg = 1) 
& 

Fig. [T] shows that robust entanglement, up to about 6cbits 
of logarithmic negativity, between two ions can be created in 
such a setup. A complete analysis of entanglement and infor- 
mation propagation through chains of ions will be detailed in 
J2I]. Here, let us just point out that such degrees of entangle- 
ment are by far out of experimental reach for quantum optics 
(where, to the best of our knowledge, ~ 1.6 is the maxi- 



mum value so far reported after state reconstruction 02511 ). The 
coupling to the bath of 7 = 10~ 6 Hz (best value considered in 
the plot) is realistic in view of the recently observed heating 
rates in ion traps I26I 271. The plot shows that this estimate is 
extremely encouraging, especially if compared to the state of 
the art for quantum optical cavities and resonators, where loss 



rates still significantly limit performances. Multipartite entan- 
glement, as well as entanglement between non-neighboring 
ions, can also be created with global control as, when the fre- 
quency is switched from u>i to u>f, all the n ions in the trap 
start interacting with each other. For instance, for three ions, 
u>i = 20 MHz and uif = 2 MHz, the initial completely sep- 
arable state evolves into a 'fully inseparable' Gaussian state 
(inseparable under any bipartition of the modes); let us denote 
the CM of this state, after an evolution time t = 5^ t _1 = 5/is, 
by <t 3 . 

Nonlocality test - Such multipartite entanglement can be put 
to use to test quantum nonlocality with massive particles. Tak- 
ing advantage of the possibility of performing parity mea- 
surements (in a single shot) and displacements, we will anal- 
yse the violation of the Bell-Klyshko inequalities lEitl on the 
three ions state with CM er 3 , by the displaced parity test 
as introduced in 12911 . In this instance, the family of (non- 
Gaussian) local, bounded, dichotomic observables is given 
by Uj{xj,pj) = b j {x j ,p j )^{-l) ft iD j {xj,pj), where D j 
and hj are the displacement and number of phonons oper- 
ators of ion j. The three observers, pertaining to the three 
ions, randomly apply two different displacements [Dj(xj,Pj) 
and Dj (x'j ,p'j)] on their ions and then measure parity locally. 
The expectation value of the operator H(R) = Ui(xi,pi) <g> 
^(#2,^2) <S> n3(x3,p3) is proportional to the Wigner func- 
tion W(R) of the composite system evaluated in the point 
R = {x u x 2 ,x 3 ,p u p 2 ,p 3 ) r : (U(R)) = (2/ir) 3 W(R) & 
Such a function is immediately determined for the Gaussian 
state under consideration: 



W(R) 



(3) 



reads: _B 3 

2>E3,i>l,P2>i J 3) 



The Bell-Klyshko inequality finally 

^ s \W(xx,X2,x's,Pi,P2,p's) + W(x 1 ,x 2 

W(x' 1 ,X 2 ,X 3 ,p' 1 ,P2,P3) - Wix'^x'frX^p'^PfrP^l < 2. 

Quantum mechanics allows for 2 < B 3 < 4. As epitomised 
by Fig. |2] regions in the space of displacements where the 
violation of the inequality is substantial and remarkably stable 
can be found. Therefore, this preliminary study opens up 
very promising perspectives concerning the violation of Bell 
inequalities with massive degrees of freedom, which would 
be a major, not yet prob ed, testing ground for fundamental 
quantum mechanics 13111 . A detailed analysis of the impact 
of imperfections and noise on such nonlocality tests will be 
presented in J^t]. 

Conclusions - We have demonstrated how the local control 
of the trapping frequencies would allow one to reproduce any 
linear optical manipulation on radial modes of trapped ions. 
We also indicated that phonon detection and homodyne detec- 
tion as well as the implementation of non-Gaussian operations 
is possible in this setting. Next, we have emphasized that, 
even restricting to global control, such manipulations enjoy a 
high efficiency in entanglement generation and low decoher- 
ence rates, along with the possibility of implementing number 
and parity measurements with current technology. The ex- 
perimental pursuit of the programme outlined in this paper 
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FIG. 2: Function B3 for pi = P2 = Ps = p[ — p'3 = 0, x[ — 
X3 = —x' 2 = — 3nm, p'2 = lnm, Onm < xi < 3nm and — 2nm < 
22 < lnm for the Gaussian state with CM cr 3 defined in the main 
text. Dimensions were reintroduced assuming v t = 1MHz and Ca 
ions. The inequality is violated in the whole displayed region, where 
the function reaches a maximum ~ 2.32. 

thus holds considerable promise, concerning both technolog- 
ical developments, such as the storage and manipulation of 
quantum information, and fundamental physical aspects, as in 
the nonlocality test for massive degrees of freedom here dis- 
cussed. 
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